Topological Invariants of Metals and Related Physical Effects 
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The total reciprocal space magnetic flux threading through a closed Fermi surface is a topological 
invariant for a three-dimensional metal. For a Weyl metal, the invariant is non-zero for each of its 
Fermi surfaces. We show that such an invariant can be related to magneto-valley-transport effect, 
in which an external magnetic field can induce a valley current. We further show that a strain field 
can drive an electric current, and the effect is dictated by a second class Chern invariant. These 
connections open the pathway to observe the hidden topological invariants in metallic systems. 
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There have been great successes in classifying insula- 
tors using topological invariants, such as the TKNN num- 
ber for quantum Hall insulators[l , and the Z 2 index for 
insulators with time-reversal symmetry [2J. The classifi- 
cation can even be extended to the more general systems 
such as superconductors (3]. In all of these systems, there 
exist gaps in their excitation spectrums, which provide 
topological protection, and resulting in robust physical 
effects. 

For metallic systems, there are also topological invari- 
ants. An apparent one is the total volume enclosed by 
Fermi surface in a non-magnetic metal, according to Lut- 
tinger's theorem. One may find other topological invari- 
ants, and in some cases, even relate them to observable 
physical effects |H[S]. For instance, in a two dimensional 
(2D) metallic system with time reversal symmetry, an 
electron moving around the Fermi cycle in the recipro- 
cal space will acquire a total Berry phase of either or tt 
per cycle. The non-trivial tt Berry phase is responsible for 
the unusual anti-weak localization effect observed in com- 
pounds such as graphene or the surface states of three- 
dimensional (3D) topological insulators [BJ. It can even 
modify the effective electron-electron interaction, giving 
rise to unconventional superconducting pairing symme- 
tries [Ml]. 

For a 3D metal, an apparent topological invariant is the 
total reciprocal space magnetic flux threading through a 
closed Fermi surface: 

<t> FS = f n-ds F , (i) 

Jfs 

where f2 (k) is the Berry curvature of the Bloch electron 
in band forming the Fermi surface at quasi-momentum k, 
defined by fi (fc) = i(V k u(k)\ X |V fc u(fc)), and \u(k)) 
is the corresponding periodic part of Bloch wave func- 
tion |101 111] . $i?5 is quantized as 

$fs = 2™ , (2) 

with n being an integer. A nonzero quantized number n 
indicates the topological nature of the metal. The topo- 
logical number n cannot be changed without changing 



the geometric topology of the Fermi surface. It is easy to 
generalize the classification to the more general metal- 
lic systems with multiple Fermi surfaces. In this case, 
one can assign a topological number to each of the Fermi 
surfaces, and the numbers cannot be changed without 
merging/splitting of the Fermi surfaces. An example of 
such topologically non-trivial 3D metals is the 3D Weyl 
(semi-)metal, in which the Fermi surfaces enclose Weyl 
nodes that give rise to the quantized reciprocal magnetic 
fluxes [T2Tfl"7] . An immediate question is how the topo- 
logical invariant can be related to observable physical ef- 
fects. 

In this Letter, we will show that the invariant can be 
related to the magneto-valley-transport property of the 
3D topological metals. We show that a magnetic field 
will induce an electric current along the direction of the 
magnetic field in each of the Fermi pockets. For a 3D 
metal that always has an even number of Dirac points 
(due to the Fermi on doubling theorem), such a response 
will generate a valley-current response, with the valleys 
corresponding to the Dirac fermions of opposite helicities. 
On the other hand, using carefully engineered strain field, 
one could create an effective magnetic field that has the 
opposite directions in the two different valleys, generat- 
ing a real electric current response. We will show that 
such response is related to a second-class Chern number 
defined on the Fermi surfaces and the real space. 

To see the response of the system to an applied mag- 
netic field, we use the semiclassical wave packet dynam- 
ics |10l fTT] . The equations of motion for the wave-packet 
center position r and momentum k in the presence of a 
magnetic field B can be written as, 

r = v n (k)-kx n n (k) (3) 
hk = -er x B (4) 

where f2„ (fc) is the Berry curvature of Bloch states of 
the n-th band, and v n (k) = de n (k) jdhk is the group 
velocity of Bloch electrons and e n (k) is the band dis- 
persion. We note that although the set of equations is 
semi-classical, it is versatile in deriving results valid in 
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Figure 1: A schematic plot to show that the electric current 
response to a magnetic field for a real 3D metal is always zero 
because the contribution from the two Dirac points always 
cancel each other. Note that for each individual Dirac point, 
the contributions from the lower and the upper bands have 
the same sign, because the normal directions (directions of 
the group velocities) of the iso-energy surfaces are opposite. 

quantum regime. Actually one can establish that the re- 
sult derived from the semiclassical equations are in gen- 
eral valid as long as the external fields are weak in the 
length scales of the lattice constants. 
It is straightforward to get: 

D n (fc) f = v n (fe) + | 0„ • fl n ) B (5) 

where D n (k) = 1 + (e/fi) (B ■ tt n ) stands for the Berry 
curvature's correction to the density of states in phase 
space We obtain a response of the electric current 
j B in the direction parallel to the applied magnetic field 
B: 



3 b = ~eJ2 / D n (k)rf(e n (*)) [*] = a B B (6) 
where the magneto-current response coefficient is defined 

by 

n J 

where <&i = f. ( fc ) =e dS ■ fl n stands for the flux thread- 
ing through the iso-energy surface at e for n-th band, 
with the normal direction of the dS being defined as the 
direction of the group velocity v n (k), and /(e) is the 
Fermi distribution function, [dk] = dfc/(27r) 3 . 

However, for a real 3D metal, the total electric cur- 
rent response is always zero. This is because the Dirac 
points, which are the source of the reciprocal space mag- 
netic field fluxes, always appear in pairs in the reciprocal 
space, and in the opposite helicities, according to the 
fermion doubling theorem |19|. As shown in Fig. [I] the 
reciprocal space magnetic fluxes from such a pair of Dirac 
points always cancel each other at all the energies. As 



a result, the magneto-current response coefficient a B is 
always zero for a real 3D metal. 

On the other hand, one can construct a non-zero re- 
sponse of valley current for systems with nontrivial topo- 
logical numbers of Fermi surfaces. To be specific, we con- 
sider a 3D Weyl (semi-)metal |I2fH7] , The simplest form 
of such a metal comprises of a pair of Weyl nodes that 
have the reciprocal magnetic monopole charge ±27rn. 
When the two Weyl nodes are well separated in the mo- 
mentum space and the Fermi level is close to the energy 
of the Weyl nodes, the electrons in the two Fermi pock- 
ets centering around the Weyl nodes can be considered 
as two independent fermion species, and we can intro- 
duce a new "valley" degree of freedom to label them. An 
external magnetic field will induced the opposite electric 
currents in the two valleys at the zero temperature: 

e 2 

j± = ±n-^ (e F - eo) B , (9) 

where ± denotes the valley index, tp is the Fermi energy 
and €q is the energy of the Weyl nodes. The correspond- 
ing valley current is: 

3v = 3+ 3 , (10) 
2e 2 

= n-j-j- (e P ~ eo) B . (II) 

Although the introduction of the valley degree of free- 
dom is based upon the separation of momentum space, 
and looks artificial, it can actually be as real as the other 
degrees of freedom such as spin, within the certain en- 
ergy and momentum scales. In particular, when n = 1, 
a pair of Weyl nodes can be mapped to the right-handed 
and left-handed neutrinos respectively, albeit in a much 
lower energy scale. When the system is clean enough or 
only has the spatially smooth disorders, the scattering 
between the two valleys is negligible, and the electron 
can maintain its valley identity for a sufficient long time. 
In this case, one expects a valley accumulation in the 
boundary of the system when there is a bulk valley cur- 
rent, as shown in Fig. [2j The detailed profile of such 
an accumulation depends on the condition of the sample 
boundary and the valley relaxation. 

Earlier attempts to utilize the valley degree of freedom 
are mainly focused on the graphene-based systems |20f - 
I23| . More recently, the monolayer molybdenum disul- 
phide (M0S2) is proposed to be an ideal material for val- 
leytronics |24l |2"5] . Our analysis shows that the Weyl 
metals can also be a candidate for the valley tronics, and 
the valley polarization can be achieved by applying a 
magnetic field. 

Next, we consider the possibility for generating a real 
electric current in a 3D metal by applying a strain field. 
A carefully engineered strain field could induce artificial 
magnetic field with its magnitude and direction depend- 
ing on the momentum |26[ I27| . Such an artificial mag- 
netic field can have the opposite directions in the different 
valleys, as observed in graphene [28 . If this can be done 
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spatial position r. Specifically, it is defined as, 




Figure 2: The magnetic field induced valley polarization. 
The magnetic field B (black arrow) will induce counter- 
propagating currents of Weyl electrons of the opposite he- 
licities (red and blue arrows). It results in an accumulation 
of right-handed (+) and left-handed (-) Weyl electrons at the 
two ends of samples, respectively. 



in a Weyl metal, from Eq. (|9]), it will generate a real 
electric current. 

The proper framework to consider the effect is to use 
the complete form of semiclassical equation: 



r = v(k, r) - *\t kr ■ r - < \f kk ■ k, 



k 



n r * • k + n rr ■ f, 



(12) 
(13) 



where v(k,r) = 
— 21m (du(r, k)/dk a \du(r, k)/dkp 



n 



kk 



de(k,r)/d(hk), ^ 



(~)fcr 

"a/3 — 

-2lm(du(r,k)/dk a \du(r,k)/dr p ), and tt rfe , 
are defined similarly [10]. Note that in the presence of 
spatial inhomogeneity due to the strain field, the Bloch 
wave function and the energy dispersion are in general 
functions of both the wave packet central momentum 
k and the central position r. As a result, the new 
components of the Berry curvatures arise. Here, ft rr 
assumes the role of that of the magnetic field B in 
Eq. Q, and can be considered as an artificial magnetic 
field. 

In most of the cases the strain field is weak, and the 
spatial gradient of the Bloch wave function is a small 
quantity. One can expand the solution of Eq. ( 12 13 I in 



the orders of the spatial gradients. To the second order 
of the spatial gradients, we find that the total electric 
current density can be written as, 



j(r) = V r x m(r) --I [dk]f(i)n kk x V r e 



(14) 



where m(r) = —(e/H) J[dk]ft kk g(e) with g(e) = 
Jde/(e), and we define a vector (fl kk ) a = (l/2)e a/ g 7 0|^. 

3>e,r is the second-class Chern-flux through the momen- 
tum space iso-energy surface S £jr with e(k,r) = e at the 
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-e aPl I dk [n^nj* + n&fi# + n r p k n k ;} , 

(15) 

where we have defined k\ and fc 2 as the two general- 
ized coordinates for parametrizing the two dimensional 
iso-energy surface S e . ri with the normal direction of the 
surface defined by v. It is easy to observe that the in- 



tegrand in Eq. ( 15 I is exactly the second Chern-form of 



Berry curvature defined in the manifold (r 1; r 2 , k2, fe), 
where r\ and r 2 are the spatial coordinates in the plane 
perpendicular to the direction of the current. 

The electric current density induced by a strain field 
has three parts of contributions, corresponding to each of 
the three terms in Eq. ( 14 1 , respectively: the first term is 



the magnetization current due to the spatial inhomogene- 
ity; the second term is the anomalous Hall current 
driven by the electro-elastic potential — V r e; and most 
importantly, there arises a topological contribution pro- 

(2) 

portional to ffr^- It is easy to see the total current across 
a cross-section of the sample I = J dS ■ j has a topolog- 
ical contribution dictated by the second Chern number 
because J dS ■ = (2n) 2 n with n being an integer, 
and we have: 



^topo. — — ~h^ ep ~ e °^' 



(16) 



We discuss the implications of Eq. ( 14 1 and Eq. ( 16 I for 
different kinds of strain fields. For a strain field created 
by the usual lattice deformation, one expects that the 
artificial magnetic field varies spatially, and there cannot 
be big spatial area with nearly uniform artificial magnetic 
field, as it costs huge elastic energy to sustain the required 
strain field |26l [27] . In this case, one cannot expect a net 
electric current induced by the strain field. This does 
not exclude the possibility of existence of a local current 
density which may in turn modify the lattice dynamics. 

A particularly interesting case is when the system has 
topological defect, e.g., a disclination, one expects an ar- 
tificial magnetic field distributing around the disclination 
line |26| . A non-zero second Chern number in Eq. (16 1 



implies the existence of a chiral conducting channel along 
the line of disclination, just like the chiral edge state of a 
quantum Hall insulator. It is important to note that the 
chiral conducting channel is in a metal , with its robust- 
ness protected by the structure of the topological defect 
and the Fermi surface topology. The construction of a 
model system to realize the novel possibility of creating 
topologically protected chiral modes in a metal is left for 
future investigations. 

In summary, we discuss the topological invariants in 
three dimensional metals. We find that non-zero re- 
ciprocal space magnetic fluxes threading through the 
Fermi surface will give rise to the magnetic-field-to- 
valley-current response, although no electric current re- 
sponse is possible due to the Fermion doubling theorem. 
We further observe that a strain field could induce an 
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electric current dictated by a second Chern number. It 
implies the possibility of existence of chiral conducting 
channel in a metallic system with appropriate structural 
topology and momentum space topology. 
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